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We report on the observation of half-quantum vortices (HQVs) in the easy-plane polar phase 
of an antiferromagnetic spinor Bose-Einstein condensate. Using in situ magnetization-sensitive 
imaging, we observe that pairs of HQVs with opposite core magnetization are generated when singly 
charged quantum vortices are injected into the condensate. The dynamics of HQV pair formation 
is characterized by measuring the temporal evolutions of the pair separation distance and the core 
magnetization, which reveals the short-range nature of the repulsive interactions between the HQVs. 
We find that spin fluctuations arising from thermal population of axial magnon excitations do not 
significantly affect the HQV pair formation dynamics. Our results demonstrate the instability of a 
singly charged vortex in the antiferromagnetic spinor condensate. 
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In a scalar superfluid, the supercurrent circulation 
around quantum vortices is quantized in units of h/m 
due to U{1) gauge symmetry [1], where h is the Planck 
constant and m is particle mass. However, when a super¬ 
fluid possesses an internal spin degree of freedom, there is 
an intriguing possibility for the superfluid to host quan¬ 
tum vortices of a fractional circulation of h/m. The su¬ 
perfluid phase interwinds with the spin orientation and a 
new relation is imposed on the supercurrent circulation 
in connection with spin texture [2]. Fractional quantum 
vortices are of particular interest in two-dimensional (2D) 
superfluidity. In the absence of long-range order in two 
dimensions [3], the superfluid phase transition in 2D is 
associated with vortex-antivortex pairing as described by 
the Berezinskii-Kosterlitz-Thouless (BKT) theory [HIS]. 
Hence fractional quantum vortices, introduced as new 
point defects, represent an interesting opportunity to ex¬ 
plore for exotic superfluid phases, possibly beyond the 
BKT physics. 

Quantum vortices having h/2m circulation, so-called 
half-quantum vortices (HQVs) have been experimentally 
observed in spinor superfluid systems such as exciton- 
polariton condensates [MSI and triplet superconduc¬ 
tors |9]. In previous cold atom experiments, HQV states 
were created with an optical method in two-component 
Bose-Einstein condensates (BECs) [10], where the two 
components are not symmetric in terms of interactions. 
Recently, a spin-1 BEG with antiferromagnetic interac¬ 
tions has been considered with great interest because 
HQVs are topologically allowed in the polar phase of the 
system HTHU]. Theoretical studies predicted an anoma¬ 
lous superfluid density jump at the phase transition in 
two dimensions [HHE] as well as a new superfluid state 
that has completely broken spin ordering [IIllIH]. 

In this Letter, we report on the observation of HQVs in 
the easy-plane polar phase of an antiferromagnetic spinor 
BEG of ^^Na atoms. Using magnetization-sensitive imag¬ 
ing, we observe that pairs of HQVs with opposite core 
magnetization are generated when singly charged quan¬ 


tum vortices are injected into the condensate. The tem¬ 
poral evolutions of the pair separation distance and the 
core magnetization reveal the short-range repulsive inter¬ 
actions between the HQVs. Dissociation dynamics of a 
singly charged vortex was previously observed in exciton- 
polariton condensates [7] but its dynamics was driven by 
spatially inhomogeneous spin-dependent potentials in the 
system. We emphasize that our system is spin symmetric 
and defect free. Thus, the observation of HQV pair for¬ 
mation clearly demonstrates the intrinsic instability of a 
singly charged quantum vortex in the spinor condensate. 

The spin-dependent part of the mean-field energy func¬ 
tional for a spin-1 spinor condensate is given as 

E, = '^-^{¥)^+p{F,)+q{Fl), ( 1 ) 

where C 2 is the spin-dependent interaction coefficient, n 
is the atomic number density, F = {Fx^Fy,Fz) is the 
single-particle spin operator, and p and q are the linear 
and quadratic Zeeman fields, respectively nans]. Here 
the external magnetic field defines the z direction. For 
antiferromagnetic interactions (c 2 > 0) and zero total 
magnetization (p = 0), the ground state of the system is 
a polar state with (F) = 0 [T9| [20]. This is a | mi? = 0) 
state along a certain quantization axis which we denote 
by a unit vector d. Depending on the sign of g, the 
condensate shows two distinctive phases: at g > 0, the 
easy-axis polar phase with fixed d \\ z, giving {F^) = 0, 
and at g < 0, the easy-plane polar phase with d F z 
and (F|) = 1. In the following, we refer to these two 
phases as polar (P) and antiferromagnetic (AF) phases, 
respectively. 

The order parameter of the AF phase can be expressed 
with a three-component spinor as 



where 'ipi is the mz = I spin component along the z 
direction (/ = 0,±1), 0 is the superfluid phase, and 
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FIG. 1: (color online). Schematic illustration of the half¬ 
quantum vortices (HQVs) in the antiferromagnetic (AF) 
spinor condensate. The superfluid phase 0 and the spin ori¬ 
entation d rotate by tt around vortex cores having nonzero 
magnetization Mz . The order parameter of the condensate is 
invariant under the operation of ^ ^ ^ + 7r and d ^ —d and it 
is continuous over the disclinations indicated by dashed lines. 


0 is the spin orientation in the xy plane, i.e. d = 
(cos 0, sin0, 0). The order parameter manifold is Maf = 
[u{i) X s^]/z2 mills], where U{1) is the gauge symme¬ 
try, comes from the rotational symmetry of the spin, 
and Z 2 arises from the invariance under the operation of 
0 — y 0 7T and (j) — y 0 T tt . 

When the windings of 0 and 0 around a quantum vor¬ 
tex are and qs in units of 27r, respectively, the single¬ 
valuedness of the order parameter requires qn ± qs to be 
integer. Therefore, quantum vortices with a half-integer 
super cur rent winding number can exist with the aid 
of spin winding. The spatial structures of the four fun¬ 
damental HQ Vs with |g^| = | are described in 

Fig. 1. When g^ + g^ = 0 (gn - qs = 0), the rriz = -1 
{rriz = 1) component has no vorticity and fills up the 
HQV core. Although a ferromagnetic core is costly for 
the AF spin interactions, the core filling is energetically 
favored as it would reduce the kinetic energy by suppress¬ 
ing the density of the circulating spin component in the 
core region. 

A singly charged vortex with {qn^Qs) = can 

be regarded as a sum of two HQVs : {qn^Qs) = 
and (±|, — |). It was predicted in mean-field calculations 
that in the AF phase a singly charged vortex state is ener¬ 
getically unstable to decay into two HQVs [2TI423] . Note 
that the two HQVs have opposite core magnetization. 
This means that if the disintegration of a singly charged 
vortex occurs in an AF spinor condensate, it would be 
identified with formation of a pair of ferromagnetic de¬ 
fects having opposite magnetization. 

Our experiment starts by generating a BEG of ^^Na 
atoms in the \F = l^rriF = 0) hyperfine spin state in 
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FIG. 2: (color online). Dissociation of a singly charged vor¬ 
tex into two HQVs. (a) Singly charged vortices are generated 
in a condensate in the polar (P) phase and then the con¬ 
densate is transmuted to the AF phase (see the text), (b) 
Optical density (OD) image of a condensate containing singly 
charged vortices, (c) Magnetization-sensitive phase-contrast 
image of the AF condensate at th = 1.5 s. The images were 
taken after 24-ms expansion by releasing the trapping poten¬ 
tial. Schematic descriptions of (d) a singly charged vortex 
state with (gn,gs) = (1?0) sind (e) a state having a pair of 
HQVs with (gn,gs) = (|, |) and (|, -|). 


an optical dipole trap [24]. The trapping frequencies 
are = (4.2,5.3,480) Hz and a typical 

condensate containing about 3.5 x 10^ atoms has the 
Thomas-Fermi radii (Rx^Ry^Rz) ~ (185,150,1.6) ym. 
For peak atom density, the spin healing length is ^s = 
h/-^2mc2n « 4.5 ym [25], which is larger than the sam¬ 
ple thickness Rz and thus, the spin dynamics in the con¬ 
densate is effectively 2D. The external magnetic field 
Bz = 30 mG, giving q/h = 0.24 Hz, and the residual 
field gradient is compensated to be less than 40 /iG/cm. 

We inject quantum vortices into the condensate by stir¬ 
ring the center region of the condensate with a repulsive 
laser beam for 10 ms, as described in Ref. [27]. Because 
the condensate is in the P phase with U{1) symmetry 
and the multiply charged vortices are unstable to de¬ 
cay [28|, it is ensured that the generated vortices are 
singly charged. The average vortex number was about 
six [Fig. 2(b)]. 

The condensate is transmuted into the AF phase by 
tuning the quadratic Zeeman field g with the microwave 
dressing technique [29H33] . First, we apply a radio¬ 
frequency pulse of 65 ys to rotate the spin from Pz to 
Tx, forming a superposition state of the ruz = ±1 com¬ 
ponents, and then we immediately turn on a microwave 
field with frequency detuned by —300 kHz with respect 
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FIG. 3: (color online). Formation of a HQV pair with oppo¬ 
site core magnetization, (a) In situ magnetization distribu¬ 
tion Mz{x,y) of the condensate at th = 1.5 s. Images of the 
HQV pairs at (b) th = 0.7 s and (c) 1.1 s, obtained by av¬ 
eraging over ten images cropped from several magnetization 
images like (a), (d) Magnetization profiles of the HQV pairs 
along the separation direction at various hold times from 
averaged images like (b) and (c). Evolutions of (e) the pair 
separation distance s and (f) the core magnetization differ¬ 
ence AMz. 
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FIG. 4: Spin fluctuations in the AF spinor condensate. In 
situ magnetization images of the condensate, containing no 
vortices, at (a) = 0 s, (b) 1.5 s, and (c) 3 s. (d)-(f) Images 

taken after 24 ms time-of-flight, where spin fluctuations are 
enhanced due to self-interference during the expansion Hi 
1^ . (g) Variance of in situ magnetization (SM^) as a function 
of the hold time th- The background noise level is subtracted. 


to the \F = l.rriF = 0) ^ |F = 2, mi? = 0) transi¬ 
tion, resulting in q/h = —10 Hz. We confirmed that the 
rriz = 0 component is absent in the microwave dressing 
using Stern-Gerlach spin separation measurements. 

The spatial distribution of the condensate magnetiza¬ 
tion is measured with a spin-dependent phase-contrast 
imaging method j34H36] . The probe light is circularly 
polarized and the frequency is detuned by —20 MHz from 
the 35 'i/2|F = 1) ^ 3 P 3 / 2 IF' = 2) transition [37]. Be¬ 
cause the phase shifts of the probe beam for the rriz = ±1 
components are opposite, the optical signal in the phase- 
contrast imaging is proportional to the density difference 
between the spin components, i.e., the magnetization 
Mz = n{Fz) of the sample. 

When the condensate containing vortices is prepared 
in the AF phase, we observe that the vortex core visibil¬ 
ity decreases in the optical density images like Fig. 2(b) 
when pairs of point defects with opposite magnetization 
emerge in the condensate [Fig. 2(c)] [38]. The magnetic 
point defects are the HQV pairs that result from the 
dissociation of the singly charged vortices. We see that 
the separation directions of the HQV pairs are random, 
which clearly indicates that the dissociation dynamics is 
not driven by external perturbation such as the residual 
magnetic field gradient but by the intrinsic instability of 
a singly charged vortex state. 

In in situ measurements of the magnetization distri¬ 
bution, HQV pairs become discernible after a hold time 


t/i ^ 0.7 s in the microwave dressing (Fig. 3). The sepa¬ 
ration distance of the pair and the magnitude of the core 
magnetization gradually increase, and the growth ceases 
after ty ^ 1.3 s when the separation distance reaches 
about So = 17.6 jum. This saturation behavior seems to 
be consistent with the short-range repulsive interactions 
between HQ Vs with opposite core magnetization [2ll|22|. 
Taking into account the imaging resolution of ~ 4 /im, 
we estimate the FWHM of a fully developed, magnetized 
HQV core to be ~ 8.4 jam that corresponds to ^ l-T^g, 
where is the average value of the spin healing lengths 
at the positions of the HQV pairs. This is in good quan¬ 
titative agreement with mean-field predictions [211123]. 

We observe that spatial fluctuations of magnetization 
develop in the condensate when it is transferred to the 
AF phase (Fig. [^. The AF condensate has two gap¬ 
less excitation modes: phonons and axial magnons. Spin 
fluctuations arise mainly from thermal population of the 
magnon mode [HITIIIT]. In our experiment, when the 
condensate is initially prepared in the AF phase by the 
spin rotation, all the spins are aligned to the Fx direction 
and this spin texture corresponds to the zero tempera¬ 
ture for magnon excitations. Thus, thermal relaxation 
of the spin temperature would subsequently occur and 
lead to spin fluctuations as observed. Because of the 
2D character of spin dynamics, we may expect further 
enhancement of spin fluctuations in our system. Spin 
fluctuations show a steady increase after an initial, rela¬ 
tively rapid growth [Fig. 4(g)], which we attribute to the 
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FIG. 5: (color online). Effect of thermal magnons on the pair 
dissociation, (a) Vortices are generated in the AF condensate 
after a dwell time tv . Thermal spin fluctuations increase with 
tv [Fig. 4(g)]. (b) Pair separation distance s as a function of 
the hold time th after vortex generation. The data in Fig. 3(e) 
are displayed together and labeled as = 0 s. The solid lines 
are the lines of so(l — with so = 17.6 /xm, fit 

to the data: (to,Ts) = (0.03,0.29) s for = 1 s (blue) and 
(to,Ts) = (0.44,0.35) s for = 0 s (red). 

heating effect of the microwave dressing [33] . 

The dissociation of a singly charged vortex involves 
spin texture formation as well as magnetized core devel¬ 
opment, as depicted in Fig. 2(e). Therefore, magnons 
that are spin wave excitations might affect the pair dis¬ 
sociation dynamics. To investigate the possible effect of 
thermal magnons, we measure the temporal evolution 
of the pair separation distance at various spin temper¬ 
atures of the condensate. To vary the spin tempera¬ 
ture, we exploit the aforementioned thermal relaxation of 
magnons and let the condensate dwell in the microwave 
dressing for a time ty before generating vortices by stir¬ 
ring [Fig. 5(a)]. Here we assume that HQ Vs cannot be 
directly generated by the stirring because the moving op¬ 
tical potential induces only density perturbations in the 
condensate. 

We find that the increasing rate of the pair separa¬ 
tion distance is not significantly affected by thermal spin 
fluctuations. This implies that once a singly charged 
vortex is split, the subsequent dynamics of HQV pair 
formation is mainly driven by the repulsive interac¬ 
tions of the HQ Vs. The characteristic time scale of the 
pair formation is estimated by fitting a growth model 
so(l — to the experimental data, which re¬ 

sults in Tg ^ 0.3 s, where td is the hold time after vortex 
generation [Fig. 5(b)]. We note that the value of is 
roughly comparable to the time scale set by the spin- 
dependent interaction energy C 2 n/h ^ 10 Hz. The model 
fitting to the data for > 0.5 s and the data of the pre¬ 
vious experiment [red circles in Fig. 5(b)] gives to ^ 0 s 


and 0.5 s, respectively. This seems to suggest that ther¬ 
mal spin fluctuations at ty > 0.5 s are large enough to 
quickly split the vortex to a certain separation distance 
beyond which the repulsion between the HQ Vs becomes 
prominent. 

In conclusion, we have observed HQ Vs with mag¬ 
netized cores and confirmed the instability of a singly 
charged vortex in the AF spinor condensate. We will 
extend this work into a 2D regime [42j [43], where a 
pair superfluid state without spin ordering is predicted 
to exist at finite temperatures HZllll]. The spin order¬ 
ing of the system might be probed spin-sensitive Bragg 
spectroscopy [44l [45] or matter wave interference meth¬ 
ods [46]. 
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